Observation of an ultracold plasma instability 
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We present the first observation of an instability in an expanding ultracold plasma. We observe 
periodic emission of electrons from an ultracold plasma in weak, crossed magnetic and electric fields, 
and a strongly perturbed electron density distribution in electron time-of-flight projection images. 
We identify this instability as a high-frequency electron drift instability due to the coupling between 
the electron drift wave and electron cyclotron harmonic, which has large wavenumbers corresponding 
to wavelengths close to the electron gyroradius. 
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Ultracold plasmas (UCPs), produced by photoioniz- 
ing a sample of laser-cooled and trapped atoms, have 
extended neutral plasma parameters by over two orders 
of magnitude, to electron temperatures below 1 K [l^. 
Studies of UCPs have primarily concentrated on temper- 
ature measurements 0, IH, and expansion dynamics 
05 H 0] , and recent work has identified a stable collec- 
tive mode 0. A signature of the collective and nonlin- 
ear nature of plasmas is the existence of plasma insta- 
bilities, perturbations that grow exponentially to large 
amplitudes and dominate plasma dynamics. Much of 
the quest for fusion energy involves control and suppres- 
sion of plasma instabilities [9]. This universal dynamics 
occurs in all kinds of plasmas, including space plasmas 
[iqI . [ill . [12I. dusty plasmas [13], magnetically confined 
plasmas [14, and Hall thruster plasmas for spacecraft 
propulsion [li. El, Ell, El • 

In this work, we present the first observation of a 
plasma instability in an expanding UCP. By applying a 
small magnetic field (~ 2 G) perpendicular to an applied 
electric field (~20 mV/cm), we observe periodic pulsed 
emission of electrons from an expanding UCP, with a fre- 
quency range from 50 to 500 KHz. Using a time-of-flight 
electron projection imaging technique [19], we image the 
electron spatial distribution by extracting them with a 
high- voltage pulse and accelerating them onto a position- 
sensitive detector. We observe that electron projection 
images split into two or three lobes in the ExB direction, 
coincident with the observation of periodic electron emis- 
sion signals. This provides strong evidence for a plasma 
instability in the expanding UCP due to the electrons 
drifting relative to the ions across the magnetic field. A 



high-frequency electron drift instability [17|, ll8| quanti- 
tatively matches our observation, which has a frequency 
lower than the electron gyrofrequency and a short wave- 
length on the order of the electron gyroradius, due to the 
coupling between the electron drift wave and a harmonic 
of the electron cyclotron frequency. 

Our production of UCPs is similar to our previous work 
[H, which we briefly summarize. We cool and trap about 
a few million metastable Xenon atoms in a magneto- 



optical trap (MOT). The neutral atom cloud has a tem- 
perature of about 20 /iK, peak density of about 2 x 10^^ 
cm~^, and Gaussian spatial density distribution with an 
rms radius of about 300 /im. We then produce the plasma 
by a two-photon excitation process, ionizing up to 30% 
of the MOT population. One photon for this process is 
from the cooling laser at 882 nm, and the other is from a 
pulsed dye laser at 514 nm (10-ns pulse). We control the 
ionization fraction with the intensity of the photoioniza- 
tion laser, while the initial electron temperature Tgo (typ- 
ically around 1-500 K) is controlled by tuning the 514-nm 
photon energy with respect to the ionization limit. 

The ionized cloud rapidly loses a few percent of the 
electrons, resulting in a slightly attractive potential for 
the remaining electrons, and quickly reaches a quasineu- 
tral plasma state. It then expands with an asymptotic 
velocity (50-100 m/s) caused by outward electron pres- 
sure [5], and maintains a roughly Gaussian density pro- 
file during the lifetime of the UCP [20]. There are two 
grids about 1.5 cm above and below the plasma to pro- 
vide a small electric field E (~5-50 mV/cm) so that elec- 
trons leaving the plasma are guided to a microchannel 
plate detector. We can apply a longitudinal magnetic 
field 5|| parallel to E, or a transverse magnetic field 
B± perpendicular to E. The applied magnetic fields are 
turned on before the two-photon excitation process. The 
black curve in Fig. [T] is the typical electron emission sig- 
nal from a freely expanding UCP for Tgo = 100 K. The 
signal consists of a prompt peak (initially escaped elec- 
trons) followed by a region of little electron loss where 
the quasineutral plasma state forms. This is followed by 
a long ~150 /is loss of electrons, interpreted as the de- 
cay of the plasma as electrons spill out of the potential 
well, which gets shallower as the plasma expands. The 
electron signal with a small longitudinal B\\ (~ 12 G) 
looks similar to that without a magnetic field except for 
small enhancement and some changes in the expansion 
dynamics (expansion dynamics in a large B\\ is studied 
by time-of-flight projection imaging technique 0). 

The electron emission signal with a small transverse 
magnetic field B± is dramatically different from that with 
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FIG. 2: The emission frequency as a function of E/B±; each 
symbol for varying B± from 1 to 3.2 G and constant E; differ- 
ent symbol for varying E from 7 to 33 mV/cm and constant 
B± . The emission frequency is linearly dependent on the elec- 
tron drift velocity Vd. 



FIG. 1: Electron emission signals for different E and B±. 
(a)-(c) constant E and different B±; (d)-(f) constant B± and 
different E. The red curves (with large periodic emission) are 
single shot; the blue curve (with much less periodic emission) 
is the average of 40 shots, the black curve (with large prompt 
peak) is the electron signal without a magnetic field. 



a small B\\ . We observe periodic pulsed electron emission 
even with a very small B±^ as shown in Fig.[T^-[T]3. As we 
increase B± from zero, the electron signal starts to have 
oscillations at about 0.8 G. The periodic emission appears 
at about 30-50 fis after the formation of the plasma with 
frequency of several hundred KHz. The three traces in 
each panel of Fig. [1] correspond to individual single real- 
izations of the UCP. Note that the emissions have similiar 
character for each shot, although the phases are random. 
As we continue to increase Bj_ to about 2-2.5 G (Fig.[T]3), 
the amplitude of the emissions gets larger, comparable to 
the prompt peak in the absence of a magnetic field. The 
frequency decreases to 50-100 KHz and the prompt peak 
gets even smaller. The background electron emission sig- 
nal almost vanishes except for few peaks at about 3-3.5 
G (Fig.Hfc). 

The frequency and amplitude of the periodic elec- 
tron signal also depends on the applied electric field E 
(Fig. [TJi-lTf). The frequency increases and the amplitude 
of the emissions decreases as we increase E, which is sim- 
ilar to the case of decreasing B±. Using fast Fourier 
transform (FFT) of the electron signals in Fig. [H we can 
extract the emission frequency as a function of E/B± 
(electron drift velocity Vd), as shown in Fig.[2l The emis- 
sion frequency depends linearly on the drift velocity with 



a slope of 10.2 ±0.2 m-^. 
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FIG. 3: Electron projection images in weak crossed magnetic 
field and electric field. All the units of the images here are in 
pixels, and one pixel is about 150 /xm. 

By using a time-of- flight projection imaging technique 
[191] . we image the electron spatial distribution onto a 
position-sensitive detector (a phosphor screen attached 
to a multichannel plate detector), and record images with 
a CCD camera. Figure [3] shows the electron projection 
images at different delay times after the formation of the 
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UCP for different B± and E. All the images are single 
shot, with little variation from shot to shot. In the ab- 
sence of a magnetic field (Fig. [3^), the electron images 
show a Gaussian density profile. Note that the electron 
image size slowly decreases during the lifetime of the 
UCP because of the strong Coulomb force of the ion cloud 
on the electrons, electron loss and electron Coulomb ex- 
plosion effects, unlike ion images (whose sizes decrease in 
the first 20 /is due to the strong Coulumb explosion ef- 
fect of the dense ion cloud, and then increase afterward, 
reflecting the real size of the plasma) As we in- 

crease B± to about 1 G (Fig. [3)3), the electron images 
start to split into two lobes at about 30 /is, coincident 
with the observation of large periodic electron emission 
signals. At about 2 G, we observe up to 3 lobes in the 
electron images between 30 /is and 70 jas (Fig. [3^1 and[3]i). 
The extra electron lobes are in the ExB± direction, and 
they show up in the other side of the main electron cloud 
if we change the sign of B±. We did not observe any 
changes in the ion images and ion current in the cross- 
field configuration compared to those without a magnetic 
field. This is quite surprising, as one typically expects the 
electron distribution to closely follow the ion distribution 
for a cold neutral plasma. In order to confirm that the 
extra lobes of the electron images are coming from the 
plasma itself, rather than dynamics during the flight to 
the detector of the projection imaging technique, we took 
a series of images at different extraction parameters such 
as the high-voltage amplitude and width and the accel- 
erating voltages on the middle and front grids (located 
between the plasma and the detector, and strongly affect 
the sizes of electron images). We always observe the sim- 
iliar results as Fig. [3l implying this is a good measure of 
the electron distribution in the plasma. Using a simple 
model of the electron space charge effects during transit 
to the detector, we find the separation between the main 
cloud and the extra lobes before time of flight to be about 
0.5-1(7^, where at is the plasma size at delay time t. That 
is, the extra lobes are located inside the plasma. 

There are a multitude of different instabilities in plas- 
mas, but given our parameters we can constrain possible 
choices to a small number. We identify the periodic elec- 
tron emission signals as well as the extra lobes in the 
electron projection images as a signature of a plasma in- 
stability due to electrons drifting relative to ions across 
the magnetic field. The crossed magnetic and electric 
fields will drive the electrons to drift with a velocity Vd 
(Yd = E/B±) in the E x B± direction. The ions are 
unmagnetized, not affected by the small Bj_ (several G), 
due to their large mass (the ion gyroradius is much larger 
than the UCP size and the ion gyro-period is much longer 
than the UCP lifetime). In our UCP the electrons are 
unmagnetized in the first ~30 /is because the electron 
gyroradius is about the same order of magnitude as the 
UCP size and the electron collision rate is higher than 
the electron cyclotron frequency. As the UCP expands. 



the plasma size gets larger, the plasma density gets lower 
and the electron temperature gets smaller (due to various 
cooling mechanisms) Q. The electrons become magne- 
tized at about 30-50 fis after the formation of the UCP 
because the electron collsion rate starts to be less than 
the electron gyrofrequency and the plasma size at that 
time is about a factor of 10 larger than the electron gyro- 
radius, which is consistent with our observation that the 
periodic emissions start at the same time. The frequency 
is in the range from 50 KHz up to a few hundred KHz, 
which is much larger than the ion cyclotron frequency 
fci (~12 Hz at 1 G) and much less than the electron cy- 
clotron frequency fee (~2.8 MHz at 1 G) and the electron 
plasma wave fpe which is about 10-20 MHz at the delay 
time of about 30-50 /is. The periodic emission signal is 
roughly independent of plasma density (or plasma fre- 
quency) as shown in Fig. [T] (note that the plasma density 
drops by a factor of 8 from 50 /is to 100 /is). There are 
periodic emissions with different frequencies (Fig. [TJd and 
[TJi), which also indicates that we have plasma instability 
(mode switching) in the cross-field configuration. 

A candidate instability is the high-frequency electron 
drift instability, which has been studied in Hall thrusters 
(a type of plasma-based propulsion systems for space- 
crafts) experimentally [l^ and theoretically 13, [l^ 
in the past few years to explain the transport of electrons 
across the magnetic field lines. A 2D fully kinetic model 
of the Hall thruser developed in ref. [13] has demonstated 
that the large drift velocity at the exhaust of the thruster 
was the source of an instability that gives rise to plasma 
turbulence and could induce a significant current across 
the magnetic field. It is a high-frequency electron drift 
instability with frequency fr lower than the electron cy- 
clotron frequency fee and short wavelength close to the 
electron gyroradius Tee, which is studied in ref. [18|. 

The theory of high-frequency electron drift instabil- 
ity is developed from the dispersion equation of electro- 
static waves in a hot magnetized electron beam drifting 
across a magnetic field with unmagnetized cold ions, and 
is closely related to the ion- acoustic- wave instabilities in 
other plasma systems, such as non-specular radar me- 
teor trail s [Tol , the ionosphere of the earth Tl|, [l2[ , dusty 
plasmas [is] and magnetic pulses [1?], but they usually 
restrict the analysis to the cases where the drift veloc- 
ity Vd is much smaller than the electron thermal veloc- 
ity Veth' For high-frequency electron drift instabilities in 
Hall thrusters and here, Vd is about the same order of 
magnitude as the Veth^ and much larger than the mag- 
netic field gradient drift velocity and the density gradient 
drift velocity. 

Assuming an electrostatic field perturbation ^ = x 
exp(ik • r — icot) and a Maxwellian electron distribution 
function with mean velocity Vd and temperature Tg, the 
dispersion equation can be written in the limit kz = 
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FIG. 4: Numerical solutions of the ID dispersion relation 
of equation (1) as in [18], but for our condition, (a) repre- 
sents the envelopes of the real (frequency) and the imaginary 
(growth rate) as a function of kyVd/oOpe- (b) the correspond- 
ing ones for small wavenumbers. The electron cyclotron fre- 
quency is equal to 0.2oUpe, and Veth/Vd = 0.5. 



about 9-14 m~^, in agreement with the measured slope 
of Fig. [2] (10.2 m~^). The maximum of the growth rate 
is reached for kyVd/ujpe = 1.2, and the corresponding 
wavelength is about 0.5 mm, close to the electron gyro- 
radius (about 0.2 mm for B± = 2 G and Tg = 2.7 K). 
We can also see the transitions from stability to instabil- 
ity whenever kyVd/oOpe is close to a cyclotron harmonic 
nQce/^pe- The growth rate reaches a maximum and then 
decreases sharply between each cyclotron harmonic, and 
is separated by stable regions. The frequency is several 
orders of magnitude below the growth rate except in the 
vicinity of the maximum. 

The frequency roughly linearly depends on the drift 
velocity (below kyVd/ujpe = 1), but is independent of the 
plasma frequency (i.e. plasma density) before the fre- 
quency reaches the maximum value (both axis are scaled 
by electron plasma frequency in Fig. [4^), which explains 
the lack of density dependence seen in Fig. [TJ We mea- 
sure the growth rate by suddenly applying the electric 
field at different delay times, and find the periodic emis- 
sion occurs within 1-2 jas after application of the field, 
corresponding to a growth rate of ~ 100-150 KHz for 
B± = 2 G and £^ = 20 mV/cm, which is about the same 
as the frequency 



(where kz is the wavenumber along B±) [it 
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xy is the wavenumber in the plane perpendicular 
to B±] A|) = kBTe/mcOp^ is the Debye length; Qce is the 
electron cyclotron frequecy; cOpe is the electron plasma 
frequency; ky is the wavenumber in the E x B± direction, 
and b = kyV^^^/ft'^. The functions In are modified Bessel 
functions of order n. 

In order to find the dependence of uOr and coi on the 
wavenumber ky^ we have solved numerically equation (1) 
by taking a further approximation to neglect the pertur- 
bations perpendicular to the drift motion {k^ = 0, the 
wavenumber along E) because the electron drift velocity 
Vd is about an order of magnitude larger than the density 
gradient drift velocity. Figure [4] shows the frequency and 
the corresponding growth rate of the unstable modes as a 
function of ky. Figure [4^ and[4)3 are the envelopes of the 
solutions and the first 5 modes, respectively. The values 
of fr for the maximum growth rate correspond to fre- 
quencies ranging from 90 KHz to 140 KHz for fpe = 21 
MHz, B± = 2 G, E = 20 mV/cm and Te = 2.7 K at 
delay time of about 30 /is (where Tg is from Q), which 
agrees with the measured frequencies at = 10^ m/s in 
Fig. m The ratio of the frequency to the drift velocity is 



Although we see large scale density structures in the 
electron images, we do not have a model of the micro- 
scopic electron motion that produces the periodic pulses. 
Note that the applied magnetic field should suppress elec- 
tron detection, since it is transverse to the detector di- 
rection. This is seen in the strong suppression of the 
prompt electron emission peak. Nevertheless, we detect 
large pulses of electrons, presumably due to large elec- 
tron trajectories that extend past the grid (1.5 cm below 
the plasma) where the large acceleration field can direct 
the electrons to the detector, even in the presence of the 
transverse magnetic field. 

In summary, we have observed large periodic electron 
emission and splitting of the electron distribution into 
two or three lobes from an expanding UCP in the pres- 
ence of crossed magnetic and electric fields. We iden- 
tify them as a signature of high-frequency electron drift 
instability due to the electrons drifting relative to the 
ions across the magnetic field. We calculate the unstable 
mode frequencies and growth rates by solving ID dis- 
persion equation and find good agreement. The large 
scale changes in the electron spatial distribution remain 
a mystery, as does the exact mechanism that leads to the 
emission of electrons. This work shows that UCPs will 
continue to provide an interesting place to study funda- 
mental plasma physics phenomena. 

This work was supported by the National Science 
Foundation PHY-0714381. 
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